Information propagation through quantum chains with fluctuating disorder 
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We investigate the propagation of information through one-dimensional quantum chains in fluctuating 
external fields. We find that information propagation is suppressed, but in a quite different way compared 
to the situation with static disorder. We study two settings: (i) a general model where an unobservable 
fluctuating field acts as a source of decoherence; (ii) the XX model with both observable and unobservable 
fluctuating fields. In the first setting we establish a noise threshold below which information can propagate 
ballistically and above which information is localised. In the second setting we find localisation for all 
levels of unobservable noise, whilst an observable field can yield diffusive propagation of information. 



Quantum lattice models as frequently encountered in the 
condensed-matter context or in quantum optics obey a kind 
of locality: once locally excited the excitation will travel 
through the lattice at a finite velocity. For spin models 
this speed at which information can propagate is gener- 
ally limited by the Lieb-Robinson bound [1] which says 
that there is an effective "light cone " (or "sound cone ") 
for correlations, with exponentially decaying tails, whose 
radius grows linearly with time The importance of un- 
derstanding information propagation in interacting quan- 
tum systems was only understood recently when it was 
exploited to establish the Lieb-Schultz-Mattis theorem in 
higher dimensions [3]. In generalising the proof of this 
breakthrough result an intimate link between the speed 
of information propagation and the efficient simulation of 
these systems has been revealed fllH]. 

The argument underlying the Lieb-Robinson bound re- 
lies only on the ultra-violet cutoff imposed by the lattice 
structure and is therefore very general. Hence there are 
some situations where the Lieb-Robinson bound is not the 
best available: when we know more about the structure of 
the interactions it should be possible to construct tighter 
bounds — an intuition which has been borne out in Refs. 
10,0] where it was shown that for the XX model with static 
disorder there exists an effective "light cone" whose ra- 
dius grows logarithmically with time [6]. In this example 
Anderson localisation IH] for the tight-binding model was 
combined with a Lieb-Robinson type argument to supply 
the stronger bound on the multiparticle dynamics. 

Interacting quantum lattice systems are also of great in- 
terest in quantum information science as they could be used 
as quantum wires to join together different parts of quan- 
tum computers or quantum devices [9, 10]. When used in 
this manner quantum spin chains become quantum chan- 
nels and the standard Lieb-Robinson bound then supplies 
a fundamental lower bound on the time required to trans- 
mit information: the receiver must wait for a time at least 
proportional to the length of the chain before the message 
arrives. For the disordered XX model mentioned above the 
bound of Ref. |@] shows that the receiver must wait for a 
time which is exponential in the chain length before infor- 



site 



FIG. 1: Schematic illustration of the different "light cones" for 
the different regimes: (a) the (lower) linear one for ordered sys- 
tems; (b) the (upper) logarithmic one for static disorder; (c) the 
(middle) diffusive one for dynamic disorder. In each regime in- 
formation is strongly attenuated outside the associated light cone. 



mation about the message can be extracted, thus rendering 
such systems essentially useless as quantum channels. 

In this work we obtain bounds on the propagation of 
information through interacting quantum spin chains in a 
fluctuating disordered field, a setting which is ubiquitous in 
the condensed matter context [11]. It will become appar- 
ent that — in some situations — information propagates dif- 
fusively; in other situations the system becomes localised 
in the sense that information can propagate effectively by 
at most a constant number of sites. 

We study two models: (i) a spin chain with general 
nearest-neighbour interactions in a noisy field fluctuating 
both in strength and direction; (ii) an XX spin chain in 
a field fluctuating in strength only. To solve the dynam- 
ics of both these models we find a master equation which 
describes the time evolution. We then explicitly calculate 
an improved Lieb-Robinson bound for the first model be- 
fore finding various correlation functions for the second 
model. This allows us to identify several different possi- 
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ble regimes of information propagation: either localised or 
ballistic propagation for the first model and either localised 
or diffusive propagation for the second. 

Master equation — We derive a master equation for each 
of our models with time-dependent Hamiltonian 



(1) 



where a^jCr^ , are the Pauli matrices, is the intrin- 
sic system Hamiltonian and the second term describes the 
interaction with the noisy field. For the first model a runs 
over the three directions x,y,z whilst for the second model 
the field direction is fixed (a = z); j runs over all lat- 
tice sites and ^f{t) are the field strengths at time t, which 
must be distributed according to a probability distribution 
with finite second moment (e.g. a Gaussian distribution, in 
which case ^" = dB"" describes a Brownian motion Bj). 

After time t an initial state p(0) evolves to p{S^, t) = 
U{C,t)p{0)U{^,ty where the propagator U{l,t) is 
given by the time-ordered exponential U t) = 

Texp (— i H{^, s)ds). We now describe the derivation 
of the master equation for the ensemble averaged density 
operator p{t) = E^p(^, t) where represents an average 
over all possible realisations of the disorder. We split the 
time evolution up into AI small time steps and at each step 
we couple each site to ancillas or meters — one for each di- 
rection x,y,z (only one meter per site is used for the model 
with fixed field direction). When each meter is initialised 
in a Gaussian state, one can take the small time step limit 
M — > oo, following the steps of continuous-measurement 
theory in |12[|. resulting in the master equation 



d,p{t) = -i[p{t),Ho] - jY.Mt),a] 



(2) 



where 7 = cE(^")^, c > 0, is a constant measuring the 
strength of the disorder. As before, for the second model 
we omit the sum over a and simply set a = z. The 
above derivation is valid when all times scales of interest 
are longer than the typical fluctuation timescale. 

General spin chain model — In this section we derive a 
Lieb-Robinson bound for a general spin chain of length n 
in an unobservable field whose strength and direction fluc- 
tuate independently on each site. The Hamiltonian is given 
by Eq. ([Xl with intrinsic Hamiltonian equal to a sum of gen- 
eral nearest-neighbor interactions Hq = hj , where hj 
acts on sites j and j + 1. Note that this model is general 
enough to include all the standard quantum lattice systems, 
including the ferromagnetic and antiferromagnetic Heisen- 
berg models and also frustrated systems. 

A Lieb-Robinson bound governs the speed of "light" or 
"sound" in the lattice model. It is commonly expressed as 
an upper bound on the Lieb-Robinson commutator. 



CB{x,t) 



sup ■ 

Ax 



[A.,Bm 



where is an operator acting non-trivially only on site 
X, B{t) = E^W{C,t)BU{C,t) is an ensemble aver- 
age over all possible realisations of the operator B in the 
Heisenberg picture, and || • || is the operator norm. While 
the Lieb-Robinson commutator is not directly observable, 
any bound for it easily translates to a bound on all ob- 
servable time-dependent two-point correlation functions 
{AB{t)) - {A){B{t)). Hence, the Lieb-Robinson bound 
is a convenient bound to express the decay of correlations. 

In the Heisenberg picture B obeys an equation simi- 
lar to the master equation derived in the previous section, 
namely dtB{t) = i[Ho,B{t)] - lYlj'Dj{B{t)) where 
the decoherence operators are defined as 'Dj{B{t)) = 
Y.o\^'^M-.B{t)\]. Defining = E,Pj(^)and 

using a Taylor expansion we arrive at 

\\[A,.B{t)]\\ < {l-^e^n)\\[A.,,B{t) + ie[H^,B 

+ e-i\\[A^,SnB{t)-V{B{t))]\\+0{e^ 

We deal with each of the non-negligible terms on the right 
hand side in turn. Making repeated use of Taylor expan- 
sions and the unitary equivalence of operator norm it is 
easy to see that the first term can bounded as || [A^, + 
ie[H,,Bm\\ < e\\[[Ho,A,],B{t)]\\ + \\[A,,Bm + 
O(e^). In order to bound the second term, we note the fol- 
lowing two identities: (i) D/j(cr") = Sdj^kO'k all a 7^ 



and Vk{a^ 



0, where cj° 
0, 
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= I; (ii) for all /? / 

i€{0,x,y,z} 



(3) 



If we define B^ = cr"^ (X) • • • a"" to be a tensor product 
operator indexed by the vector a = (oi, . . . , a„), it is 
straightforward to use the above identities to show that 

\\[A,,8nB^-V{B^)]\\ 

<Ek^x.,e{o,x,y.z}\U^^'^kBc.a^k]\\ (4) 
= 8(n-l)||[A„Sj|| 

We can extend this result to general B by linearity: 
B{t) = '^^Ca{t)Ba. This implies that the sec- 
ond term obeys the bound (|4l with i?„ replaced by 
B{t). Combining the above equations allows us to 
bound the time derivative dfCsix, t) < —8'jCb{x, t) + 
sup^ II [[i?o, Aa;], ||. It remains to rewrite the 

off-diagonal term, the right-most term containing the dou- 
ble commutator. Note that [i/o) ^£c] = 2||i/o|| H^^H^ 



x + 1) 



with V = Ea,/3e{0.x.v,.}«-l^: 

where lu^'^^l < ||y|| < 1. By defining the vector CB(t) = 
(CB(l,i), . . . ,CB{n,t)) and employing the above expan- 
sion for V we are able to rewrite the bound on the time- 
derivative of the Lieb-Robinson commutator as 

dtCeit) < (-8(7 - 8\\Ho\\)I + 32\\Ho\\R)CB{t) 

where the inequality is to be understood componentwise 
and Rj^k = ^j.k+i + ^j+i,k- This has solution 



CB{x,t) < e 



-8(7-8||Ho||)t 



„32||Holl-Rt 



.C7s(j,0).(5) 
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FIG. 2: Illustration of the dynamics of the correlation function 
Cj,i{^,t) = {^l\aj{0)a\{t)\Q) generated by Eq. © for a 50 site 
chain for a specific realisation of the fluctuating field £,j{t). The 
horizontal axis is site number j and the vertical axis is the time 
(for 7 = 0.05, 7 = 0.1, 7 = 0.2 and 7 = 0.5). 



This is the main result of this section: a new Lieb-Robinson 
bound for our class of models. The key idea here is that 
the first term can (provided the noise 7 is large enough) 
exponentially suppress information propagation, whilst the 
second term can cause ballistic propagation of information. 
We now analyse the interplay between these two effects: 

Discussion of the new bound — Noting that ||i?|| = 2, 
it is clear that if 7 < 16||//o|| then our bound is expo- 
nentially growing and the original Lieb-Robinson bound 
is superior to ours: it limits us only to ballistic propaga- 
tion of information. If, however, 7 > 16j|i^o|| then the 
right-hand-side of our new bound (|5]) is negligible if ei- 
ther (i) \j — k\ > tK^ for some constant > or (ii) 
t > t, = log{C/e)/r where C = Ej <^b(j,0) and 
r = 87 + 128||i?o||- Equivalently, the right hand side of 
the new bound is non-negligible only when |j — A;| < ^ete, 
which is to say that non-negligible correlations can propa- 
gate by at most a constant number of sites. 

XX spin chain — In this section we abandon the Lieb- 
Robinson commutator in favor of various correlation func- 
tions as they are easier to calculate in this specific scenario. 
In the previous section we studied a noisy field fluctuating 
in both strength and direction but here we consider a fluc- 
tuating field oriented in the z-direction. For convenience 
we study the XX model, although one can expect similar 
results to hold for other spin chain models with nearest- 
neighbor interactions and a fixed field direction. We begin 
by applying the Jordan Wigner transformation to map our 



system of qubits into a system of spinless fermions 

HiC, t)=y] (a]a,+i + a]a, 1 + ait)a]aA (6) 



The sum is over an infinite chain or a finite ring (in which 
case we impose periodic boundary conditions). In order 
to study the propagation of information through the sys- 
tem we introduce disorder dependent correlation functions 
which are the probability amplitude for a fermion to hop 
from one site to another in a given time interval 

c,,,{^,t) = {n\a^U{(,t)al\n) 

where \Q) is the vacuum. We further define ensemble 
averaged correlation functions Cj,k{t) = E^Cj^fc(^,t) = 
{n\aj{t)al\n) where a^{t) = E^W {(„t)a^U{C,t) is 
the ensemble averaged annihilation operator in the Heisen- 
berg picture. Using Eq. (|2l we form a differential equation 
fortr[aj(t)/)], yielding 5ta^(t) = -z(a^^i(t)+aj_i(t))- 

7a^(t). Introducing a (t) = {ai{t),a2{t), . . .) allows us to 
rewrite this as a vector differential equation with solution 
a{t) = e~^*e*^*a(0). This in turn implies that the ensem- 
ble averaged correlation functions are 
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(e 



)j,k 



An analysis similar to that performed previously for the 
Lieb-Robinson bound reveals that all ensemble averaged 
correlation functions are exponentially small excepting 
those for which |j — A;| < n{e,j) = const. In other words, 
the information can — on average — propagate by at most a 
constant number of sites. In parallel with the results of the 
first model we have localisation of information but in con- 
trast we have no noise threshold. See Fig.|2l 

It must be noted that this result holds only when we av- 
erage over the disorder: specific realisations of the disorder 
could lead to less localised dynamics. In order to see how 
far the dynamics of a specific realisation of the noise can 
stray from the averaged dynamics found above we restrict 
ourselves to the infinite chain. We introduce the discrete 
squared position operator and the momentum operator 
00 

(where \ j) denotes the state vector of a single excitation at 
site j). The Heisenberg picture time derivatives obey equa- 
tions similar to Eq. e.g., dtX^{t) = i[HQ, x^{t)] — 
-f£{x'it)) where £{M)/2 = M - Y., 
eliminates the diagonal entries of M. When the initial state 
is a single particle at site 0, /9(0) = |0) (0|, then 



2t e-27t _ 1 
— + ^ 



--■■fi7,t) 



Recalling the trivial bound Pj,j{t) < 1, using translational 
invariance, and noting that {x^){t) = ^jppj.jit) allows 
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us to conclude that when we place the particle initially at 
site k, p(0) = \k){k\, then 

<min{l,/(7,t)/|j-/c|2} (7) 

One can define ensemble variance correlation functions by 
v,,k{t) = E4|c,-fe(C,t)p - |E^c,-fe(e,t)p. A little alge- 
bra verifies that E^jc^ t)p are precisely the diagonal 
elements Pjj{t) when p(0) = \k){k\; consequently the 
Vj^k{t) also obey the bound (|7]). 

This bound implies that it is unlikely that we will stray 
far from the averaged dynamics given by Cj,k{t) for small 
times (2t/7)^/^ <C \j — k\. By employing Chebyshev's 
inequality which states that for a random variable X with 
mean /x and finite variance and for any positive real 
number k then P {\X — /u| > na) < 1/ k?, we can con- 
clude that there is an effective light cone whose radius 
grows proportionally to /(7, t)^^'^ (2t/7)^/'^. In other 
words, information can — on average — propagate by a dis- 
tance proportional to the square -root of the time elapsed, 
reminiscent of a classical random walk lll3n . 

To aid understanding of the above result, we express 
the averaged motion in the Heisenberg picture to find that 
wave packets stop in their motion: since [Hq,p\ = 
and £{p) = 2p, the Heisenberg picture momentum obeys 
= E^tr[p(t)p(0)] = e"'>'*(p)(0), so an initial exci- 
tation "localizes" in momentum space exponentially fast. 

Mixing properties of generic fluctuating disorder — 
What is the steady state the dynamics converges to on aver- 
age for very long times? We now finally link the dynamics 
under disorder to unitary designs showing that 

generically one will obtain the maximally mixed state on 
average. Consider the following master equation for a spin 
chain (with 7 > 0) 



d,p{t) = -i[p{t),H,]-^Y.^[p{t),X,lX, 



(8) 



by -i[i?o,^a] = J^fS^a.pBp, 

• • • (8 o""" as before. The opera- 



Define the matrix F 
where = cr^^ 
tors Xj can without loss of generality be taken to be of 
the form Xj = xaj + yfT|, and we write I = {a : Oj £ 
{x, y} for some j}. If the operators governing the local 
disorder have full rank and if the submatrix of F of entries 
Fa (3 for which a / and /3 € / is a matrix of maximum 
rank, then linit^oo = linit^,, E^p(e, t) = (1/2)®" 
lUql . meaning that the state becomes maximally mixed. Al- 
most all Hamiltonians H^) and local fluctuations have this 
property. Time evolution under fluctuating disorder there- 
fore approximates a unitary 1-design |14, [Tsll arbitrarily 
well for large times and the system "relaxes" on average. 

Summary and outlook — In conclusion we have studied 
various regimes of fluctuating on-site disorder (noise) in 
quantum spin chains and we have identified several regimes 
of information propagation. These have been compared to 
the propagation regimes for similar models with both static 
disorder and no disorder. We have found that in some in- 
stances (namely those with low noise levels or fixed field 



direction) that the localisation due to fluctuating disorder is 
weaker than that caused by static disorder; conversely we 
have identified other regimes (those with high levels of un- 
observable noise) for which the localisation is stronger than 
that caused by static noise. This highlights the complicated 
nature of disordered quantum systems. Indeed, such local- 
ization under time-dependent disorder should also be di- 
rectly observable in experiments with cold atoms, exploit- 
ing similar ideas of fluctuating speckle potentials as used 
in the seminal experiment of Ref. [17]. One could even 
think of using fluctuating disorder in cold atom systems as 
a "disorder filter", shaping traveling wavepackets in their 
form or letting them stop. 

It is a simple matter to extend the results above by replac- 
ing the intrinsic Hamiltonian Hq with a time dependent one 
Hf){t): the same master equation is obeyed. This would al- 
low us to implement quantum logic gates on neighboring 
qubits hence providing a means with which to accomplish 
quantum computation. The important point is that this pro- 
vides a more realistic model of noise than has been hith- 
erto used by many authors, including Refs. [18]: the usual 
scheme is to assume perfect gate operation followed by an 
instantaneous error; our model allows one to analyse sys- 
tems where the qubits are subject to noise and decoherence 
at all times and in particular whilst a gate is being applied. 
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